A MultiBaker Map for Thermodynamic Cross-Effects in Dynamical Systems 
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A consistent description of simultaneous heat and particle transport, including cross effects, and 
the associated entropy balance is given in the framework of a deterministic dynamical system. This 
is achieved by a multibaker map where, besides the phase-space density of the multibaker, a second 
field with appropriate source terms is included in order to mimic a spatial temperature distribution 
and its time evolution. Conditions are given to ensure consistency in an appropriately defined 
continuum limit with the thermodynamic entropy balance. They leave as the only free parameter of 
the model the entropy flux let directly into a surroundings. If it vanishes in the bulk, the transport 
properties of the model are described by the thermodynamic transport equations. Another choice 
leads to a uniform temperature distribution. It represents transport problems treated by means of 
a thermostatting algorithm , similar to the one considered in non-equilibrium molecular dynamics. 
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I. INTRODUCTION 

Irreversibility in transport models based on dynamical 
systems with only a few degrees of freedom have become 
a subject of intensive recent studies (l]dTl). They illus- 
trate how macroscopic transport coefficients are related 
to the properties of the microscopic dynamics. It is a 
remarkable discovery that in chaotic dynamical systems 
a rate of irreversible entropy production can be defined 
|,|§J|,|l|-|o|. This development opens the possibility 
of requireing for a consistent dynamical-system modelling 
of an irreversible process the derivation of both the trans- 
port equations and the entropy balance. In our approach 
we shall observe this constraint. 

Many models PHlO| were originally designed to rely on 
equations of motion, where transport is induced by an ex- 
ternal field and the (average) work done by the field on 
the systems is taken out by a so-called Gaussian ther- 
mostat. This approach is commonly implemented with 
periodic boundary conditions. It has extensively been 
tested numerically, but gives rise to conceptual problems 
in the interpretation of the entropy (or heat) flux, since 
it only allows to address the global entropy balance of 
a macroscopic system and has no boundaries where en- 
tropy fluxes can be let into the environment. The aim of 
the present article is to investigate a class of dynamical 
systems taylored to describe simultaneous particle and 
heat transport, driven by appropriate boundary condi- 
tions and an external field. We work out local entropy 
balance, and identify conditions under which the model 
can be consistent with non-equilibrium thermodynamics. 

Earlier low dimensional models are devoted exclusively 
to either particle (mass, charge) transport J|,|],|2] |7| or 
to heat conductivity (2^j2^]. These transport processes 
are driven by a single thermodynamic force. It is of ba- 
sic interest, however, to understand how thermodynamic 
cross-effects generated by the presence of two indepen- 



dent driving forces can be described in the framework of 
dynamical system theory (as for many particle systems 
see, e.g., [|o|). The case of thermoelectric phenomena in 
which the driving forces are (i) the temperature differ- 
ence and (ii) the electric field and/or a density gradient 
is illustrated by Fig. |l|. The cross effects imply that the 
temperature gradient contributes to the electric current, 
and the density gradient to the heat current. Recently, 
the authors of the present paper suggested an elementary 
model to study these effects j3l]]. Here, we generalize it 
to explore the conditions under which a consistency with 
thermodynamics can be found. 
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FIG. 1. Graphical illustration of the transport process con- 
sidered. A system of spatial extension L is attached to reser- 
voirs inducing particle and heat currents due to the differ- 
ences in the densities (n) and temperatures (T) at the two 
ends, as indicated by the arrows. Along the chain heat can 
be exchanged with a thermostat. 

Part of the above mentioned problems with the en- 
tropy balance has recently been clarified ||l6|-|l8|| in the 
framework of multibaker maps modelling quasi one- 
dimensional particle transport at constant temperature. 
These models are given in terms of the time evolution of 
the (multibaker) phase-space density in a correspond- 
ing two-dimensional single-particle phase space, which in 
this case consists of a band of length L. The fixed width 
represents a phase-space variable p in addition to the spa- 
tial coordinate x along the band. The multibaker map 
defines a discrete dynamics (cf. below), which is one-to- 
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one on the band, but does not necessarily preserve the 
volume locally. A recent paper by Tasaki and Gaspard 
p]| shows that analogous results can be obtained with 
area-preserving maps by making the width of the band 
position dependent. This varying height was connected 
to changes in the (potential) energy, and does not appear 
as a driving force independent of the density gradient or 
the external field. 

Multibaker maps have no natural momentum variable 
conjugated to the spatial coordinate x. Therefore, we 
characterize the thermodynamic states by two indepen- 
dent fields. Besides the phase-space density g, a new field 
w is introduced, whose dynamics describes the evolution 
of the kinetic- energy per particle,, i.e., gw corresponds to 
the kinetic-energy density plf . 

To make contact with non-equilibrium thermodynam- 
ics, the time evolution of average densities in regions of 
small spatial extension along the x-axis is considered. 
They will be called coarse-grained densities. Transport 
equations in the form of differential equations are ob- 
tained in a continuum limit, where the spatial resolution 
of coarse graining is much smaller than the linear size L 
of the system and the time unit r of the discrete dynam- 
ics is much shorter than that of the macroscopic relax- 
ations. In this macroscopic limit the phase-space density 
g and the kinetic-energy per particle w are related to 
the particle density n and to the local temperature T, 
respectively. 

It is not obvious that a deterministic dynamical sys- 
tem as simple as a multibaker can fulfill all constraints 
required for consistency with thermodynamics |3^| — not 
even when taking the macroscopic limit. The thermody- 
namic entropy balance relates the time derivative of the 
entropy density s to the entropy production a^ 1TT ' per 
unit volume and to the entropy flux $ 



of the particles, and v e i is the drift velocity due to an 
external electric field E, which is related to a e i by 



<7elE 



en 



(4) 



II is the Peltier coefficient, and a the thermoelectric 
power (or Seebeck coefficient). Thermodynamic cross ef- 
fects are manifest in (||), and the corresponding Onsager 
relations imply a relation between the transport coeffi- 
cients a and II. 

Since the entropy plays a central role in these rela- 
tions, one has to choose an appropriate entropy con- 
cept for the multibaker. One obvious candidate is the 
Gibbs entropy S^ G ' defined with respect to the phase- 
space density g. Due to the ever refining phase-space 
structures, the chaotic dynamics generates from every 
smooth initial density, this entropy never becomes time- 
independent, not even in a macroscopically steady state. 
In contrast, a corresponding entropy S, whose definition 
is based on coarse-grained densities, is thermodynami- 
cally well-behaved and its value per unit length is the 
analog of the entropy density s appearing in (Q). This 
entropy will be called the coarse-grained entropy. The 
irreversible entropy production of arbitrary steady and 
non-steady states was identified as the time derivative of 
the difference between the Gibbs and the coarse-grained 
entropy 



Although the thermodynamic relation (2b) requires 
the entropy flux to be the divergence of the entropy cur- 
rent, we allow for deviations from thermodynamics in 
that we do not exclude the presence of an additional term. 
This term is interpreted as the consequence of a thermo- 
stat, which can locally remove or release heat, leading to 
an additional entropy flux $( thcrmostat ). I n suc h cases 
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We consider thermoelectric phenomena induced by par- 
ticles of charge e in a transport process along the x axis. 
In a system which is translation invariant perpendicular 
to the x axis, the quantities causing entropy changes can 
be expressed as |33] 



r (i") 



e 2 f 



CTelT 



+ A 



d x T 
T 



(2a) 
(2b) 



Here o~ e i and A denote the electric and the heat conduc- 
tivity, respectively, and 



J = nv e i 
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W = A^ 



(3a) 
(3b) 



are the particle and the entropy current densities, re- 
spectively. In Eq. (ph, \x c denotes the chemical potential 



We call a system thermostatted, whenever $( thermostat ) 
differs from zero. Depending on the details of the model, 
we are able to study both non-thermostatted and ther- 
mostatted systems, and in the latter case we shall be 
able generate arbitrary stationary temperature profiles. 
To our knowledge, these features have not yet been ex- 
plored in non-equilibrium molecular dynamics simula- 
tions 

The paper is organized as follows. In Sect. |l|we revisit 
the thermodynamics of irreversible processes by rewriting 
the expressions of entropy production as well as of parti- 
cle and entropy currents in forms amenable to a compar- 
ison with the results of multibaker maps. Subsequently, 
(Sect. H|) we introduce the multibaker map, and discuss 
the time evolution of the phase-space density and the 
kinetic-energy density. The Gibbs and coarse-grained en- 
tropies and their dynamics are studied in Sect. [V. Subse- 
quently, in Sect. M the macroscopic limit of the obtained 
expressions is taken. Conditions on the baker dynamics 
to make it consistent with thermodynamics are explored 
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in Sect. VI. A short discussion is devoted to thermostat- 



ted cases (Sect. VII), and our main results are summa- 
rized in Sect. |VIII| . The paper is augmented by two ap- 
pendices. App. |A| is devoted to a formal definition of the 
map and the resulting time evolution of the densities. In 
App. ^ it is shown that the macroscopic results do not 
depend on the prescription for coarse graining. 



II. NON-EQUILIBRIUM THERMODYNAMICS 

In this section we recall the thermodynamic descrip- 
tion of transport induced by two independent driving 
fields [^2|. The most general situation is treated, which 
comprises the presence of an external field, as well as 
gradients in the particle density and the temperature. 



A. Thermodynamic forces and currents 

We consider a system of particles of charge e in a non- 
moving background. Due to an (electro-) chemical poten- 
tial gradient and a temperature gradient, both a particle 
and an energy current is flowing through the system. Let 
j and j^ u \ respectively, denote the density of these cur- 
rents in a frame of reference fixed to the background. In 
this setting the number density of particles n and the 
density u of the total energy are locally preserved. (The 
density u also contains the potential energy due to an 
external field and the kinetic energy of the ordered mo- 
tion.) 

In order to derive the entropy balance in a region of 
fixed volume, we start from the conservation laws 



d t n = -Vj, 



(6a) 
(6b) 



and express the time derivative of the entropy density 
per unit volume s 



8 t s = -Vj (s) +a a " 



(6c) 



in terms of the entropy-current density and the irre- 
versible entropy production er( irr ) . 

Considering a system in local equilibrium, the Gibbs 
relation 



T ds = du — /idn, 



(7) 



holds with T and /i denoting the local temperature and 
electro-chemical potential, respectively. To find the time 
derivative of s, we write the local temporal change of (M) 



t s = ^8 t u-t dtn 



(8) 



in the form of (^) by identifying 



as the total entropy current, and 

r (irr) = j(u)yl --V7 M 



T jV^ = -(/ s >VT + jV M )l (9b) 



as the irreversible entropy production. 

Equation ( pbj ) shows that the currents j and j^' 
are conjugate to the thermodynamic forces — V/i/T and 
— VT/T, respectively. Therefore, in the spirit of non- 
equilibrium thermodynamics, these currents can be ex- 
pressed in terms of the forces as 



T V/x VT 

J — -■L'll-j, ^12-jT, 



21 22 f-^-f 
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L21 

£1: 



-J, 



(10a) 
(10b) 
(10c) 



where t he la tter ex pres sion for j 1 -^ was obtained b y in - 
serting (10a) into ( 10b ). Moreover, by inserting ( |l0b| ) 
into (9b), one expresses the irreversible entropy produc- 
tion as 



(irr) 



L n L 22 - L 12 L 21 f VT\ 2 f_ 



L\ 2 — L 2 i VT . 
r tttJ- 



(11) 



After using the Onsager relation L 12 = £21 j this expres- 
sion for the irreversible entropy production is a quadratic 
form, which takes only nonnegative values provided that 
the matrix of kinetic coefficients is positive definite, i.e., 
the well-known condition L11L22 — £12^21 > is fulfilled 
fill. 



B. Identifying transport coefficients 

It is worth expressing the kinetic coeffcients Ly by 
means of directly measurable quantities. The total 
electro-chemical potential can be split as jj, = \i c + e<j>, 
where \i c is the chemical part, and (j) is the electric po- 
tential. Since E = — and j e i = ej is the electric 
current, we find that Lu is proportional to the electric 
conductivity er e ; > 0: 



L u 



(12a) 



In the absence of a particle current (i.e., for j = 0) , 
TjW provides the heat current, so that in view of (10c) 



L11L22 — L\ 2 L 2 



(12b) 



(9a) 

where A > is the heat conductivity. 



3 



At zero particle current and constant chemical po- 
tential, a temperature gradient induces an electric field, 
which is conventionally written as ceVT, where a is called 
the thermoelectric powe r (or the Seebeck coefficient). 
Consequently, from (fOa) one finds 



L12 — 



acr el 
~fe' 



(12c) 



Finally, in a sytem without temperature gradients, the 
entropy current due to the presence of an electric current 
ej amounts to Tl ej/T, where II is the Peltier coefficient. 
Hence, Eq. (|lOc|) implies 



Ifov 



>21 



(12d) 



Using the phcnomenological coefficients (|12|) we write 

(Tel 



3 



-(V/i + eaVT), 



•W - -A— — ' 



(13a) 
(13b) 



Note that the Onsager relation L 12 = £21 makes the 
Peltier and Seebeck coefficients connected as 

n = To.. (14) 

Substituting Eqs. ( |l2|) into (|ll]) one thus derives (pa]). 

C. Relating transport and diffusion coefficients 

It is worth replacing the chemical potential in the ex- 
pressions for the currents and entropy production by the 
density n and temperature T. We write 



De 2 



OeZ 



-Vn - ?VT, 



where the diffusion coefficient is defined as 



D = 



and 



dn 



dT 



(15) 



(16a) 



(16b) 



is a quantity of dimension entropy per particle. By means 
of the drift velocity v e \ [E q. (Q)] one can rewrite the par- 
ticle current density (13a) as 



where 



nD 

j = Vein - DVn - k^^VT 



nD a e i 
k— = — (ea-?). 
1 e z 



(17) 



(18) 



The entropy current then takes the form 



J 



'A- knD^pJ ^ + ^(v el n - DVn), (19) 



and in view of (jlj), (16a) and (|l§| ) one obtains 
' du c 



ell = kn 



dn 



-T 



dT 



(20) 



Note that in the entropy current ( |T9| ) the coefficients in 
front of the v e in and the —DVn terms coincide. 

We finally mention that by taking the limit E, e — > at 
constant ell, ea, and a e i/e 2 , the thermoelectric problem 
is formally mapped onto the problem of thermal diffusion 
in a binary mixture. In that case n stands for the concen- 
tration of one of the diffusing materials, and the quan- 
tity kD is the thermal diffusion coefficient |^| . Based on 
this analogy, we consider kD in ( |l8| ) as the thermal diffu- 
sion coefficient of charged particles in the thermoelectric 
problem. 
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FIG. 2. Graphical illustration of the action of the multi- 
baker map on the phase space (x,p). (a) The mapping is 
defined on a domain of N = L/a identical rectangular cells of 
size a x b, with boundary condition imposed in two additional 
cells and N+1. (b) After every time unit r each cell is 
divided into three columns, which are squeezed and stretched 
to obtain horizontal strips of width a. The average values of 
the fields pjx, p) and w(x,p) on the cells (strips) [cf. Eqs. 
( ^l| ) and (|25|)] are given on the margins. Iteration of this rule 
defines the time evolution. 
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III. THE MULTIBAKER MAP 

Baker maps are known to be prototypes of strongly 
chaotic systems |3(|. Multibaker maps are a generaliza- 
tion, where a spatially extended system is modelled by 
a chain of mutually interrelated baker maps, in order 



to model transport ||22 , g3j , g5jjlg |18| . |21| ] via the dynamics 
of the (multibaker) phase-space density g. The single- 
particle phase space modelled by the multibaker map 
consists of N identical cells of width a and height b, which 
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are labelled by the index m (Fig. ||a). After each time 
unit t every cell is divided into three columns (Fig. ||b). 
The right (left) column of width ar m ( al m ) is mapped 
onto a strip of width a and of height 1 ( bl m -\ ) m 
the right (left) neighbouring cell. The middle one pre- 
serves its area s~m = s m . The map is globally phase-space 
preserving, i.e., s m + r m + l m = + l m + = 1. A 
formal definition is given in App. [A|. The map mimics 
the time evolution of a many-particle system with weak 
interactions in a single particle phase-space |Io|. The 
equilibrium equations of state will turn out to be those 
of a classical ideal gas. 

The (x 7 p) dynamics of the multibaker map is consid- 
ered to be a microscopic dynamics in the sense that it 
is deterministic, chaotic, and mixing. It drives two fields 
which generate ever refining phase-space structures. 

Besides these fields, we also consider coarse-grained 
densities obtained by averaging over the cells. The cell 
width a is considered to be the smallest length scale 
where a thermodynamic description applies. In the spirit 
of local thermodynamic equilibrium, the local averages 
of the microscopic variables characterize the thermody- 
namic state in the cells. The temporal evolution of the 
coarse-grained versions of g and w is consequently ex- 
pected to describe an approach towards a steady state, 
where the the coarse-grained fields do no longer change 
in time (in contrast to the fully resolved fields, which 
approach closer and closer towards fractal distributions) . 
To emphasize that coarse-graining is taken over the cells, 
the coarse-grained density will also be called the cell den- 
sity. 

The dynamics of earlier multibaker models is the same 
for all cells. There can be inhomogeneities in the den- 
sities, but the evolution equations are kept translation 
invariant. Here, we relax this constraint by allowing the 
transfer rates of cell m to depend on the coarse grained 
fields in cell m and its neighbors. This dependence mim- 
ics the effect of the thermodynamic driving force due to 
for instance a local temperature gradient. It induces 
an m-dependence of the parameters. Since all calcu- 
lations can be performed without ever referring to the 
form of these dependencies, we will not yet specify them 
but start from a map with a general set of space and 
time-dependent parameters. Their form will be fixed a 
posteriori by comparison with thermodynamics. 

Due to the self similarity of the multibaker dynamics, 
the local transport and entropy balance can be worked 
out by a calculation considering one time step only. A 
discussion of more general prescriptions for the coarse 
graining is relegated to App. 0. 



A. Evolution of the phase-space density 

Thermodynamic transport equations describe the time 
evolution of the phase-space density g and the kinetic- 
energy density w. For explicit calculations of their time 



evolution we always start with the constant values g m 
and w m in every cell m. This is convenient from a tech- 
nical point of view, and does not lead to a principal 
restriction of the domain of validity of the model as it 
will be demonstrated in App. [b|. After one step of itera- 
tion the densities will be piecewise constant on the strips 
i = R,S,L defined in Fig. ||b. Due to the conservation 
of the particle number, they are 



/ _ 'm— 1 
Qm,r = — • Qm—Xi 
I'm 

6m, s = Snn 
I Im+l 



(21) 



In 



Qm+1- 



The factors r m _i/7\ n and l m +i/lm give rise to local con- 
traction or expansion of phase-space volumes. 

The coarse-grained density g' m after one time step is 
the average of the contributions (|2l]) on the different 
strips: 

g' m = (1 - r m - l m )g m + r m _ig m _i + l m+ ig m +\- (22) 

Multiplying the equation by r" 1 and introducing the cur- 
rent 



. _ oh 

Jm " yrnQm ^m+l£?m+l ) 
T 



(23) 



through the right boundary of cell m, Eq. (|22j) appears 
in the form of a continuity equation 



KQm - Qm) 



Jr, 



J, 



(24) 



It can be seen as the discrete counterpart of (|6a|). Note 
that by definition the current through the left boundary 
of cell m is the same as the current through the right 
boundary of cell m — 1. Other types of currents associ- 
ated with cell m will also be defined as a flow through 
the right boundary. 



B. Evolution of the kinetic-energy density 

The (x,p) dynamics does not imply any constraint on 
w. Its time evolution can be chosen according to phys- 
ical intuition. In contrast to the particle density, we 
consider the kinetic energy per unit volume g m w m as 
a non-conserved quantity. Besides a contribution from 
the particle flow, the new values w' m i on the strips 
i = R,S,L will therefore contain terms characterized by 
a local source strength q m , which accounts for a local 
heating: 



Qm,r ^m,r 
6m,l w 'm,l 



r m -i 



g m -i w m -i [1 + rq m ] , 



Qrn 

[1 + rq m ] , 



(25) 



In 



■ Qm+1 w m+ i [1 + rq„ 
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The source term q m is taken constant in every cell since 
more general coices only effect terms which drop out in 
the macroscopic limit. The particular form of q m will 
only be specified later, in order to demonstrate that there 
is a unique choice for q m , where the time evolution of 
the kinetic energy can become consistent with thermo- 
dynamics. Other choices for the source term q m lead 
to a non- vanishing entropy flux <j>( therinostat ) ) anc i can be 
considered to characterize a thermostat in the spirit of 
non-equilibrium molecular dynamics. 

An update of the kinetic- energy density can be calcu- 
lated similarly to the update of the particle density. In 
cell m the average value Q' m w' m after one time step is 
obtained by averaging the different contributions on the 
strips in cell m [cf. Eq. (p5|)1, yielding 

Bm w m — [QmWm{l - T m — l m ) + r m -lQ m -lW m -l 

+l m +iQ m +iW m+ i] [I + rq m }. (26) 

For a fixed set of transition probabilities r m , s m , l m and 
q m = 0, Eq. ( p6| ) amounts to a passive advection of the 
field to by the deterministic dynamics. The possible de- 
pendence of the transition probabilities on differences of 
the coarse-grained w and the presence of the source q m , 
however, make the advection nonpassive. 

Equation ( p6| ) can also be rewritten in the form of the 
discrete balance equation 



tD I 
~aJ \ 
tD 



b {Qm w m - QmWm) 



In 



J 



(qw) 



Aew) 

Jm—l 



1 + rq n 



The first term on the right-hand side characterizes the 
source strength of the field g m w m per unit time t, and the 
second one is the discrete divergence of the gu>-current 



J 



ab , 

T 



(<?») = ™f r w 



Qrn + lW m+1 ). 



(28) 



Since gw plays the role of a kinetic-energy, we consider 
j{gw) £ an ener gy current. 



C. Diffusion and drift 

The local transition probabilities r m , s m and l m gov- 
ern the evolution of the coarse-grained densities g m and 
w m . In view of the master equations (|2^) , the cell-to-cell 
dynamics of the model is equivalent to the dynamics of 
random walkers with fixed step lenght a and local tran- 
sition probabilities r m and l m over time unit r. Such 
random walks are characterized by the local drift v m 
and diffusion coefficient D: 



J m ; 



a 
2r 

r J- I = n 



Hence, the transition probabilities r r , 
expressed as 



(29a) 
(29b) 
and l m can be 



It, 



2D 



f _ av m \ 
\ 2D J 



(30a) 
(30b) 



We allow in the present paper only for a location de- 
pendence of the drift v m , but keep the diffusion coef- 
ficient spatially homogenous. The m-dependence of the 
drift should be a consequence of the inhomogeneity of the 
kinetic-energy (temperature) gradient along the chain. 

In spite of the freedom we still have in specifying v m , 
these definitions already allow us to write the currents in 
a form very close to their thermodynamic counterparts. 
The current for the phase-space density appears in the 
form 



Jm = -^{g m V m + Qm+lV-m+l) - OD — 



Q<n 



(31a) 



a discrete version of (|3^). Similarly, for the gw-current 
one obtains 



j^f) = w m j m - bg m+ iD- 



1 - 



■Um+l 

2D 



(31b) 

which describes an advection of the w fields by the parti- 
cle current and a contribution from the discrete gradient 
of the kinetic energy. 



(27) 



D. Parametrizing phase-space contraction 



Due to the condition r m +l m = r m +l m , one can express 
T\n and l m in an analogous way as (pm by introducing 



an additional parameter e via 

. tD / av 

r m = - T [l + e 



a 
tD 



1 



2D 
2D 



(32a) 
(32b) 



Using ( p2| ) and (|30|) to evaluate r m+ i — l m , one easily 
verifies that 



r m +i 



, • (33) 

is constant along the chain. This number fully character- 
izes the phase-space contraction of the multibaker map. 
In harmony with the common use in the dynamical sys- 
tems literature, we will also say that e characterizes the 
dissipation. It will be kept constant when taking the 
macroscopic limit. 

There are two special values of the parameter e which 
are worth mentioning. When e takes the value e — 1, the 
phase-space dynamics is locally area preserving (r^Xi = 

r m ,l m = lrn+i)- For e = -1 (rQTi = l m+ ul m = r m ) we 
call the resulting phase-space dynamics time reversible, 
since the initial area of any small region is recovered af- 
ter taking an arbitrary closed path along the chain (in 
points differring from the initial one, however, the area 
is in general different from the initial one). 



G 



IV. ENTROPIES AND THEIR TIME 
EVOLUTION 



denotes Boltzmann's constant). We write g*(w) in the 
form 



A. Gibbs entropy S (G) 



Q*M = 



f(w(x,p)) 



(35) 



For the generalized multibaker map, the Gibbs entropy where Q* is a constant reference phase-space density, and 
is defined in terms of the phase-space density g{x,p) and / is a dimensionless function. The actual form of f(w 



the kinetic-energy density w(x,p) as 



S^ G ' = —k B / dxdp g(x 1 p) hi — 



g(x,p) 



will be determined below by the requirement of consis- 
tency with thermodynamics. At the moment we assume 
only that it is sufficiently smooth to expand it to second 
order in its argument w. 

(G)' 

The Gibbs entropy S m after one time step can be 
where g*(w) is a reference density, which depends on w, expressed by making use of and J25| ) for the three 
and, through it, also on the phase-space coordinates (k B columns of cell m: 



g*(w(x,p)) ' 



(34) 



S<f )' = -k B ab 



(1 - r m - l m )g m In ( —f{w' m s ) 



-r m -ig m -i In 



' m— 1 Qm—X 



m,r)j +'m+lftn+ll n ^ ~ ^ f( w m,l)^j 



(36) 



After inserting the update ( p2| ) for the phase-space density, substracting Sm = —k B abg m ln(g m f(w m )/g*), and 
rearranging terms, one finds 



o(G)' _ a(G) 



(Q m - Qm) In —f(w m ) J + g m In 



— 1 Qm — 1 



in 



Tm-l gm-l f{ w 'm,r) 
T^n g-m f{w ') 



f( W ' m ,s) 
f{Wm) 



+ Im+lQm+l In 



lm+1 Qm+1 f( w 'm,l) 
IZ Qm f(w') 



,•(») _ 

Jm 7m- 1 



(j)( tho 



4>. 



(37) 



This can be interpreted as a balance equation for the Gibbs entropy. The temporal change of comprises two 
contributions: the divergence of an entropy current 



;( s ) = 



j m k B In ( ^f{w m ) J + k B a m+1 gm+l In 



gm+l /(w m+ l) 
gm f(w m ) 



and a flux into the thermostat 



hermostat) 



k B 



dm ln 



"^m-lgm-l In 



f[Wm) 

1 rn — 1 



r m /(O/K-l) 



+ ^m+lgm+1 in 



^m+1 



(38a) 



(38b) 



In general, this decomposition is not unique. It will turn 
out that (j)( thcrmostat ) contains terms, which can be com- 
bined to a divergence and hence transferred to the en- 
tropy current. This freedom can only be removed in the 
macroscopic limit, where in physically relevant situations 
the splitting is unique. In any case, the form of the cur- 
rent is close to the thermodynamic one ([H]): it contains 
a contribution proportional to the current density j and 
another term which by its dependence on the function 
f(w) characterizes the local kinetic-energy gradient. 



We identify the temporal change of the Gibbs entropy 
with the flux of the coarse-grained entropy. This is mean- 
ingful from an information theoretic point of view. Af- 
ter all, the Gibbs entropy characterizes the information 
encoded in the microscopic time evolution of a system. 
Consequently, changes of this entropy may only be due to 
an entropy current and to a coupling to the thermostat, 
i.e., to terms like those identified in Eq. (|37|) . 
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B. Coarse-grained entropy S 

The coarse-grained entropy S m of cell m is defined in 
an analogous way as the Gibbs entropy (53), but using 
now the cell density g m and the cell's kinetic energy den- 
sity w m 



-k B abg m In — f(w 



Qn 



er 



(39) 



The coarse-grained entropy of cell to after one time step 
is 



S' m = -k B abg' m In ( ^ f(w' m 
Q 



(40) 



It depends on the updated values of the coarse-grained 
quantities. 

In order to find the balance equation for the coarse- 
grained entropy, we make use of the argument of the 
previous subsection that the change in the Gibbs entropy 
may be interpreted as the macroscopic entropy flux. This 
allows us to rewrite the temporal change of the coarse- 
grained entropy as 



<?'—<? 



T 

(G)' _ 



s, 



( G )'\ fa _ q{G)\ 



(41) 



The first term on the right-hand side is the entropy flux 
A e S m /T through cell to, and the second one represents an 
irreversible entropy production Aj5 m /r. Hence, Eq. ( f4l| ) 
constitutes a discrete entropy balance in the form 



A e S m — S m ^ iS" m , (43a) 

AiS m — {Stti ~ S m ) — {S m ~ S m )• (43b) 

In the information theoretic interpretation of entropies, 

( G) 

the difference S m — S m measures the information of a 
microscopic system which cannot be resolved by a coarse- 
grained description. Hence, AiS m is the increase per 
time unit r of the information which cannot be resolved 
when characterizing the state of the system by coarse- 
grained densities. It is positive by construction, and ex- 
cept for a transient behaviour obtained for certain initial 
conditions (which will not be considered here), it can only 
increase. 

(G) 

Note that the entropy S m and the difference S m — S m 
might depend in general on details of the coarse grain- 
ing. This dependence drops out in the macroscopic limit, 
when calculating temporal changes (cf. App. H). All 
quantities appearing in the entropy balance ( f42| ) will turn 
out to be thermodynamically well-defined observables. 

As a last step, we discuss the explicit form of the the 
rate of irreversible entropy production. The initial condi- 
tion that the system is prepared with uniform densities in 

(G) 

every cell, implies S m = S m . The irreversible entropy 
change during one time step is therefore the difference 
between the Gibbs and the coarse-grained entropy taken 

' ( G)' 

after one time step AiS m — S m — S m ■ It can be split 
into two parts. An /-independent part, which comprises 
contributions due to the particle current, and another 
one, which is related to inhomogeneities in the kinetic- 
energy density, and hence in temperature. We write 



(42) 



(44) 



with 



where 



Ai5 (particlc) = kBab 



-g m In \-r m -xQ m -x In 

Qm 



-1 Qr, 



Im+lQm+l hi 



lm+1 Qm+1 



I, 



(45a) 



and 



A i: Sl hcat ) = k R ab 



,' In 



/KJ 



-XQn 



-1 m ~ — ; + Lm+lQm+l in- 



f(w'm,s) 



(45b) 



All terms appearing in (|||) have a proper physical mean- V. THE MACROSCOPIC LIMIT 
ing. The respective first ones characterize the change of 

entropy due to the time evolution of the coarse-grained A. Definition of the limit 
fields. The others amount to an entropy of mixing of re- 

gions in phase space with different phase-space or kinetic- The macroscopic limit implies that a << L> N > 1 and 

energy densities, respectively. Note that contributions T ig much smaller than typical macroscopic time scales 

from phase-space contraction appear only in A i 5'^ ar IC c> . (f or instance L 2 /D). Formally it is defined as 

a,r->0 (46) 
such that the spatial coordinate 
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x = am 



(47) 



C. Irreversible entropy production 



is finite. The phase-space density g m integrated over the 
momentum-like variable p becomes the local particle den- 
sity n(x). Since g m is independent of p, the integration 
corresponds to a multiplication with the vertical cell size 
b. As mentioned earlier, the field w is assumed to go over 
into the local temperature T(x) in the macroscopic limit. 
Thus, we have 



bg n 
w„ 



n(x), 
CT(x), 



(48a) 
(48b) 



where C is a constant of dimension one over tempera- 
ture. Moreover, the local drift, diffusion and the source 
strength q are kept finite while taking the limit 



D -> D, 

v m —* v(E,T,d x T) 
q m -> q(x) 



(49a) 
(49b) 
(49c) 



E denotes the external field. In the following we do not 
write out the x-dependence of the fields, of the drift and 
of the source term explicitly. 



B. Number and kinetic-energy density 

We first notice that under the assumption of smooth- 
ness the spatial dependence of the two fields can be ex- 
pressed as 



For the contribution (45a) of the particle current to the 
irreversible entropy pro duction we obtain (cf. the analo- 
gous calculation in JlT|Jlq] for details) 



(particle) 



ksClT 



Je 

nD' 1 



where 



e 1 
— nv = — 



nv — Dd x n 



(55a) 



(55b) 



agrees with the particle current (^TJ) up to the dissipation 
dependent term (1 — e)/2. 



The other contribution (45b) to the irreversible en- 
tropy production comprises the explicit dependence on 
the kinetic-energy field, can be evaluated in a straightfor- 
ward manner by Taylor expanding the function / and the 
logarithms to quadratic order around w m s . The terms 

linear in / (w m s ) / f(w m s ) exactly cancel. (Here f'{z) 
denotes the derivative of f(z) with respect to z). In non- 
vanishing order the macroscopic limit is therefore 



(heat) 



= nD(d x CTf 



f"(CT) 
J[CT) 



f f(CT) \' 
\ f(CT) J 



(55c) 



Note that the square bracket can also be written as the 
second derivative of In /. 



D. Entropy flux 



bg m ±i — *• n ± ad x n + —din, 

w m±1 -> CT±ad x (CT) + yflg(CT). 



(50) 



In order to calc ulate the macroscopic limit of the par- 
ticle current (31a) we use ( |30| ) and (|4^), to obtain 



j = vn — Dd x n. 



(51) 



From (pi|), the time evolution of the density n can be 
obtained as: 

d t n = -d x {nv) + Dd 2 x n. (52) 
Similarly, we have for the time evolution of riT 



[cf. Eqs. (|27j), ( |31b[ )1 

d t (nT) = -d x j^ +nTq 
with the nT-current 



J 



(uT) 



Tj - nD d x T. 



(53) 



(54) 



By expanding f(w m ±i) to linear order around w m , one 
finds for the macroscopic limit of the entropy current 

-j— = -j In 



nf(CT) 



Dd x n + nD d x (CT) 



f'(CT) 
f{CT) ' 

(56a) 



where n* = bg* is a reference particle density, which is 
constant in space and time. 

The macroscopic limit of the entropy flux into the ther- 



mostat (38b) is found to be 

^(thermostat 



= -nCTimiq 



f(CT) 



v 
D 



{1-ef 



d x (nv). (56b) 



It contains the spatial derivative d x (nv), which underpins 
our earlier statement that the splitting of the entropy flux 
into the divergence of a current and a flux going directly 
into a thermostat is not unique. It is natural to remove 
the derivative from < j ) ( thcrmostat ) j which leads to the en- 
tropy current 
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»•(•) 

Jl 



1 + ln 



n/(CT) 



nD d x (CT) 



f'(CT) 
f(CT) ' 

(57a) 



and to the flux 



(thermostat) 



= ~f(CT) 



v 
D 



{l-ef 



nv + eDd x n 



(57b) 



However, there are other e-dependcnt splittings, too. For 
instance, by adding (1 + e)nv/2 to the entropy current 
one obtains 

j« nfjCT) f'{CT) 

l— = -J ln ; Je + nD d x (CT) ——, (58a) 

k B n* fyCT) 



and 



(thermostat) 



= - nCT Wf) q 



l-ev. l + e n 

Tje H — nd x t>. 



2 D" 



(58b) 



This shows that the splitting of the total entropy flux into 
a divergence of an entropy current and a flux $( thormostat ) 
is in general not unique for arbitrary values of e, not even 
in the macroscopic limit. 



VI. CONSISTENCY WITH THERMODYNAMICS 

Having found the general expressions for the macro- 
scopic limit of the particle and the entropy flux, and of 
the irreversible entropy production, we are now in a posi- 
tion to make specific choices for the parameter e, for the 
yet undetermined functions q(x) and f(w), and for the 
functional v(E,T,d x T). 

Comparing ( |5l"| ) with the thermodynamic particle cur- 
rent (jl7|), we find that the drift v(E,T, d x T) must take 
the form 



Vei ~ k^d x T. 



(59) 



Since earlier we have not found any neccessity to fix v, 
this choice is obviously consistent with thermodyamics. 
It remains however to be seen if the other constraints can 
be fulfilled. 

The form of / can be fixed by observing that the term 
(d x T/T) 2 appears in the irreversible entropy production 
( pa| ) with the same coefficient as the —d x T/T term in the 
ent ropy curre nt (|3b| ). C omp arin g the /-dependent parts 
in (55c) and (|6aj) [or (57a) or (58a)] we find that this 
can only happen if 



z(lnf(z))" = -(In f(z))'. 



(60) 



The solution of this differential equation is a power law 



/(*) = z- 



(61) 



with 7 as a free constant parameter. A constant prefactor 
can be absorbed into the definition ( [j5| ) of g*. 

Concerning the value of e, there are several constraints, 
which all lead to the same unique choice, (i) The re- 
quirement to have the same coefficient in front of the 
Vein and the —Dd x n terms in the entropy current jl9| ) 

[cf. 



fixes the value of e to be — 1 [cf. jk given by (58a)]. 
(ii) A natural splitting of the entropy flux into the nega- 
tive divergence of the entropy current and a flux into the 
thermostat holds for e = — 1, too. In particular, we then 



have J] 8 



■0) 

32 



and <j)( thermostat ) 



= $ 



(thermostat) 
2_ 



(iii) 



The particle flux into the thermostat (57b) [or ( p8q )] has 
a well-defined meaning only when its second term con- 
tains the full particle current j. For all these reasons the 
only dynamics which leads to physically acceptible results 
corresponds to the choice e = — 1, which was connected 



to a time- reversible dissipation mechanism in Sect. [II D. 
For this choice also the pa rticl e-current dependent part 
of the entropy production (55a) contains the full square 
of j as required for consistency with the corresponding 
thermodynamic contribution (|2a|). 
With these choices 
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' B 



1+ln 



i(CT)- 



nD~f 



d x T 



(62) 



In view of (|3b|), we identify the thermal conductivity A 
and the Peltier coefficient n as 



A = ksnD^f, 
T e 



11 

— (CT)- f 



(63a) 
(63b) 



By this, the transport coefficients could be expressed by 
system parameters. It is remarkable that finite values 
were found for A and n although only the finiteness of 
v and D was assumed in the course of the macroscopic 
limit. 

Similarly, for the flux into the thermostat one obtains 



^ (thermostat) 



VJ 

— = jnq - —. 

k B D 



(64) 



It contains a term vj/D corresponding to the change of 
entropy associated with Joule's heating due to the drift 
v of the particles. In thermodynamics <j>( thcrmostat ) — o 
in the bulk, so that the source term q takes the form 



k B 



A 



(65) 



It describes the increase of the local kinetic energy due 
to dissipative heating. The heat thus deposited in the 
system will be transported to the boundaries by a heat 
current. 
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The entropy density obtained as the macroscopic limit 
of Eq. ^ ■ 



is 



= —ksn In 



— (CT)- 



(66) 



This implies that, up to an additive constant, ell/T is 
the entropy per particle s/n. 

Since the multibaker map describes a system of weakly 
interacting ]40| particles, it is natural to assume that not 
only its entropy function ([36|) but also its chemical po- 
tential corresponds to that of a classical ideal gas. We 
take 



[i c = k B T(l + j) 



T-. 

n 



(67) 



From these two equations of state 7 follows to be the 
specific heat at constant volume, measured in units k B . 

By substituting the chemical potential (|67]) into the 
thermodynamic expression of the Peltier coefficient ( p(i| ) 
and comparing it with the particular form obtained for 
n/T in (63b), one immediately sees that the coefficent k 
of (p9|) has to vanish. Hence, <r = ea from and from 
( |16h| ) one recovers the Onsager relation II = Tot. The 
fact that k turns out to be zero seems to be a special 
feature of the baker model with three strips. In this ver- 
sion temperature cannot move without an explicit parti- 
cle motion, and hence no thermal diffusion is expected. 

Next, we consider the heat current j( heat ) , i.e., the en- 
ergy current from which the potential energy of the ex- 
ternal field is excluded 



■{heat) = jrj(a) _|_ 



(68) 



Inserting the explicit form of currents j and [cf. (|5 
and @] into Eq. (Bal), this yields 



J 



(heat) 



-Xd x T + k B T 7 j. 



(69) 



This form is indeed consistent with thermodynamics, and 
also with the macroscopic limit of the nT-current (54). 
Multiplying the latter by fcs7, we recover (I 



It is worth poiting out that due to the definition (16a) 
of the diffusion coefficient and the particular form of the 
chemical potential (pTj), we find that 



a e ik B T — e 2 Dn. 



(70) 



Consequently, Einstein's relation holds in the multibaker 
map. One can thus express the electric conductivity by 
the diffusion coefficient D in all formulas. In particular, 
we find 



r (i") = A 



\ T 



nD 



(71) 



a formula which in the case of constant temperature has 
already been derived in earlier versions of the multibaker 
map @||. 



As further consequences of Einstein's relation, we men- 
tion that: (i) the electric drift v e i ([|) is proportional to 
E/T: 



Vet 



eDE 
~k^T 



(72) 



since the diffusion coefficient is assumed to be cons tant , 
(ii) comparing heat and the electric conductivities ( |63a| ) 
and (|7C|), we find 



A 



7 



(73) 



which implies that this ratio is independent of thermo- 
dynamical state variables. Thus, the Wiedermann-Franz 
law pq] proves to hold for the multibaker model, (iii) the 
elementary Drude theory p8| of metallic conduction pre- 
dicts the Seebeck coefficient a — U/T to be proportional 
to k B /e, i.e., to be independent of temperature or den- 
sity, which contr adic ts observation. Such a term is indeed 
present in Eq. ( 63b ), but its second term also predicts 
a specific state dependence. Thus, the present model 
turns out to describe certain features of transport more 
realistically than the classical Drude model, although it 
cannot be expected to give a microscopically realistic the- 
ory of thermoelectric phenomena, (which contain essen- 
tial quantum effects due to the strong degeneracy of the 
fermionic electron gas) . 

Finally, we consider the temperature equation follow- 
ing from It takes the form 



d t T = qT 



d x {\d x T ) .d x T 
k B "fn 



3- 



(74) 



which can be shown to be consistent with the general 
relation (cf. XIII. (85) of p2||) expressing the entropy's 
local time derivative as 



d t s 



d x {\d x T) 



(75) 



Here, the respective terms accounts for the heat conduc- 
tion, the Peltier and the Joul e hea ting. Substitute (66) 
and the Peltier coefficient, Eq. ( |63h| ), to see that Eq. (75) 
becomes an identity if and only if the thermodynamical 
choice of Eq. (|65|), the source term, is taken. 



VII. THERMOSTATTING 

After having identified the condition for full consis- 
tency with thermodynamics in the form of $( thcrmostat ) = 
or q = q*, we turn to a short discussion of cases where 
there can be an entropy flux into the thermostat. In 
the thermostatting algorithm of non-equilibrium molec- 
ular dynamics [Q, heat is taken out of the system in 
order to keep the temperature constant in a spatially 
homogeneous steady state, and to avoid an overheating 
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due to the permanent acceleration produced by en elec- 
tric field. In our setting this corresponds to a case with 
d x T = 0. Such a uniform temperature field is stationary 
for q — only, as follows from the temperature equa- 
tion @, so that ^(thermostat) = -k B vj/D. It is indeed 
a kind of Joule's heat, which is let into the thermostat. 
Note, however, that classical thermodynamics does not 
admit a stationary homogeneous state d x T = to be 
steady, since the temperature increases in the bulk due 
to Joule's heating. This indicates that thermostatting is 
a tool by which one can turn a preselected temperature 
profile into a steady state. After all, for every density 
profile consistent with given boundary conditions and 
a preselected fixed temperature profile T(x) there is a 
source term distribution q(x), such that the temperature 
does not change in time [cf. ([74]) with <9 t T = 0]. In all 
these cases $( thormostat ) = k B {"fnq — vj/D) is different 
from zero. This shows that the algorithm for thermostat- 
ting can be maintained even for complicated temperature 
profiles. Since the source q appears neither in the cur- 
rents, nor in the irreversible entropy production, nor in 
the constraint for the Onsager relation to hold, nor in the 
transport coefficients, the local entropy balance is consis- 
tent with every choice of q or $( th e rmostat ) provided the 
entropy flux appears in the form of (|^). In general, the 
divergence of the entropy current contributes to the 
entropy flux <E>, and the remaining part of the reversible 
change of the entropy is transferred to the thermostat. 

We conclude that although thermostatting is a devia- 
tion from classical thermodyanamics, it seems to be the 
weakest possible deviation in the sense that except for the 
form of the entropy flux it leaves all local thermodynamic 
relations invariant. It can thus be seen as an idealization 
of a physical thermostat (which can in reality only be at- 
tached to the boundaries of a system), where heat need 
not be transported spatially (to the boundaries), but can 
directly be released into the surroundings. This gives 
rise to a non-thermodynamic contribution to the entropy 
flux in the local entropy-balance equation in a general- 
ization of non-equilibrium thermodynamics [cf. Eqs. (j^) 
and (0)]. This generalization of the local balance equa- 
tion, however, does not imply at all a similarity of global 
transport properties, which in general depend on the spa- 
tial dependence of the fields. It is clear from (|7^ ) that a 
state which is steady with a given q(x) will not be steady 
with the thermodynamic choice q* (x). It will not even 
have similar density profiles. Therefore, thermodynam- 
ics and thermostatted descriptions might lead to strongly 
different results on the global level. 



VIII. DISCUSSION 

We have extended multibaker models by augmenting 
the density field g of these models by a temperature-like 
field variable w. This allowed us to address problems 
like thermoelectric cross effects requiring two indepen- 



dend thermodynamic driving fields. The model has the 
following features: 

(A) The evolution equation of w requires source terms 
reflecting the local irreversible heating of the system in 
the presence of transport. 

(B) The temperature enters the entropy through a 
kinetic-energy dependent normalization of the (phase- 
space) density. 

(C) Consistency with the thermodynamic description 
of transport is achieved for densities which are coarse 
grained in regions of small spatial extension. 

(D) Comparing the coarse-grained description with the 
microscopic one, allows us to identify all contributions to 
the local entropy balance. 

(E) The time evolution of the system can be inter- 
preted as that of weakly interacting particles, whose 
motion may only be coupled through a mean-field like 
dependence of the evolution equations on the coarse- 
grained field variables. In accordance with this, the re- 
sulting "multi-baker" gas obeys the classical ideal-gas 
equation of state. The Onsager relation, the Wiedemann- 
Fanz law and the Einstein relation can be derived, and 
expressions are found for the Peltier and Seebeck coeffi- 
cients. 

(F) The local entropy balance of non-equilibrium ther- 
modynamics can be generalized by introducing at every 
location an instantaneous flow of entropy (i.e., of heat) 
into a thermostat. When time reversibility is maintained, 
the dynamics becomes closely reminiscent to numerical 
algorithms related to Gaussian thermostats. 

(G) Dissipation and thermostatting play different roles 
in dynamical-system models for transport. The condition 
e = — 1, for time reversibility, was found independently 
of the choice of $ (thermostat) _ with tllis dissipation we 
can describe both thermostatted and non-thermostatted 
systems. 

It is remarkable that an agreement with thermody- 
namics could be achieved by this comparatively simple 
model. Indeed, strong restrictions on the choice of its pa- 
rameters were needed. The free functions (drift v, source 
strenght of heat q and normalization of densities /) had 
to be chosen appropriately, and a free dissipation-related 
parameter (e) had to be fixed to a given value leading 
to a time-reversible dissipative dynamics (even a Hamil- 
tonian, volume- preserving dynamics is excluded). With 
the given choices, however, we do not find any restric- 
tion to weak gradients. This might be a consequence of 
the baker map's piecewise linear character and of the 
Markov properties resulting from the fact that this fam- 
ily of maps admits no pruning. 

We have used a generalized concept of dynamical sys- 
tems to model open boundaries where transport can be 
induced by appropriately chosen boundary conditions. 
The time evolution of a macroscopically large number of 
independent 'particles' is considered. Consequently, not 
even in steady states the natural measure of the multi- 
baker map is relevant to calculate physical observables. 
After all, this measure is only defined if the map is closed 
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by periodic boundary conditions. Rather another mea- 
sure, the one forced on the system by the open boundary 
conditions, plays the central role. Such measures were 
first investigated by Gaspard and coworkers |j|,|22|,^3| . 

Finally, we draw attention to the fact that the present 
model differs in important features from other models of 
transport by low-dimensional dynamical systems. The 
transition probabilities (which are closely related to the 
drift and diffusion coefficients) may depend on the coarse- 
grained fields. This dependence leads to a dynamical sys- 
tem with many degress of freedom. Consequently, the full 
time evolution of the (x,p) dynamics together with the 
time dependence of the parameters can be interpreted 
as a peculiar coupled map lattice designed to closely fol- 
low transport equations. In a steady state, however, the 
parameters take time-independent values such that the 
time evolution of the particles is described by a two- 
dimensional dynamical system, a map acting on (x,p), 
which is in general lacking translation invariance. 

It is worth emphasizing that a further reduction of 
the dimension is impossible. By neglecting the p vari- 
able (i.e., when projecting the baker map to obtain 
a one-dimensional map describing the tranport of par- 
ticles along the x-direction) one finds full consistency 
with macroscopic transport equantions, but all drift- 
dependent terms disappear from the entropy balance, 
which therefore deviates from its thermodynamic form 
H ence > modelling only the transport processes via 
dynamical systems is a much easier enterprise than aim- 
ing also for a proper description of entropies. From the 
point of view of a correct entropy balance the existence 



of a phase-space variable orthogonal to the transport di- 
rection is essential. Only in this case can the fractal 
structures in the microscopic densities be followed, whose 
unresolvability leads to entropy production. It is open at 
present, however, if a dissipative dynamics is needed for 
this, since a variation of the cell size in the sense of |2l| | 
might convert contributions to the entropy production 
due to local phase-space contraction into those of mix- 
ing. 

The suggested method for modelling thermoelectric 
cross effects can be considered as a combination of a 
dynamical-system and a hydrodynamic description. Be- 
sides the appearance of a source term of the kinetic en- 
ergy, the strong mixing character of the chaotic dy- 
namics is essential, which leads to fractal phase-space 
patterns in the considered forced measure. By that it en- 
sures that irreversibility, and thus consistency with ther- 
modynamics, can be reached in a description based on a 
low-dimensional dynamical system. 
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APPENDIX A: RIGOROUS DEFINITION OF THE MAP AND OF THE TIME EVOLUTION OF 

DENSITIES 

Every point (x,p) with x 6 [0,a7V], 1 < m < N, and y € [0,6] is mapped by the multibaker map B as follows 

l m x x + a(m — 1), l m -iP for <(x — am)/a< l m 

B(x,p) = ^ s" 1 [x - al m ] + am, blrn^^p, for l m <(x-am)/a< l m + s m 

^[x - a{l m + s m )]+ a{m+\), b(l m+1 + s^Xi) + f^ip, for l m + s m <(x-am)/a< 1 

(Al) 



The dynamics of the phase-space density g(x,p;t) can 
be interpreted as the time evolution of a typical set of 
points distributed in the phase space. It is governed by 
the Frobenius-Perron equation 

g(z;t + r) = J dz' S(z - B(z'))g(z';t), (A2) 

where z denotes points (x,p) in the phase space of the 
multibaker chain. 

The boundary conditions are that g{z;t) is fixed to go 
and gN+i in the cells m — and m = N + 1, respectively. 
They are taken into account by the equation 



g{z; t + r) = g 0(N+1) / dz' 8{z - B{z')) (A3) 



for points z whose preimages lie in cell (N + 1). 

Due to the chaoticity of the map and the difference 
in the boundary conditions, the density g{z;t) becomes 
more and more irregular as time goes on (at least for any 
smooth initial distribution). Therefore, asymptotically 
the concept of density is not well defined. For t — > oo 
one assigns the measure fJ-(A) = lim^oo J. dz g(z; t) to 
any phase-space region A. 
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Similarly to the phase-space density g(z;t), the time 
evolution of the kinetic- energy density 



u> = gw, 

is described by the integral equation 



(A4) 



lu{z; t + r) = (l + Tq(z j) J dz' S(z - B(z'))cj(z' '; t). 

(A5) 

In generalization of the Frobenius-Perron equation (|A^ ), 
however, a source term q(x) is included now, which is is 
piecevise constant in the cells. The corresponding bound- 
ary conditions for w(x,p;t) are 

uj(x,p;t + T) = g ( N+ i)W 0{N+1) [l+rq(x,p)] (A6) 

for points whose preimages are in cell or N + 1, respec- 
tively 

Also in this case, w(x,p;t) is no longer well-defined 
asymptotically. Instead, the stationary kinetic-energy 
distribution should be considered as an invariant mea- 
sure v, different from /i, which assigns the weight v(A) = 
limt—.oo J A dz lu(z; t) to every region A in phase space. 

p q 



bs' T m 

mm 



bs m l m 



b L 'm+1 



R 
S 

L 

R 



L 

R 
S 
L 



R 



FIG. 3. Notations used to define the level-2 partitioning of 
a cell after two applications of the map. Every strip is labeled 
by a two-letter sequence p, q, where p,q € {L, S, R} label the 
level-2 strips, (p stands for the level-2 strips within level-1 
strips selected by the label q). To illustrate the use of these 
labels, the heights of a few strips are indicated at the left side 
of the cell. 



APPENDIX B: STRUCTURAL STABILITY OF 
MACROSCOPIC RESULTS 

In this appendix we adapt the arguments of jlTj to 
demonstrate that, after taking the macroscopic limit, 
our results are independent of the detailed prescription 
of coarse graining. One part of such a demonstration 
should be that the same results are obtained when coarse- 
graining is applied to any number m c of successive cells 



after every time step. In this case all previous results 
for quantities of unit volume remain unchanged in the 
macroscopic limit for any fixed m c , since both the length 
of the investigated region and the differences of the values 
of the coarse-grained fields become multiplied by a factor 
m c , and thus the macroscopic gradients do not change. 

We therefore concentrate on the more involved case 
of coarse graining on the cells only after every n time 
steps. From a density, which is constant in each cell, the 
dynamics generates structure on the nth level of the re- 
fining partition obtained by the nth images of the cells. 
Structural stability is demonstrated by showing that the 
total entropy production can be decomposed into a sum 
of contributions from coarse-graining on the different lev- 
els (starting with the finest structures) and that all these 
contributions are the same in the macroscopic limit. This 
is equivalent to showing that coarse graining may also to 
be applied to a finer partitioning within the cells, with- 
out affecting the resulting thermodynamic relations. For 
sake of conciseness we immediately explore the physically 
relevant case f(z) — z^ 1 [cf. Eq. (|6l|)1, suppress the nor- 
malization constant g* of the density and Boltzmann's 
constant ks, and we only work out the case of n=2. Al- 
though this corresponds to just the simplest possibility, 
it indicates the strategy to be followed when discussing 
further refinements 139(1 . 



1. Preliminaries 

We first extend our notations in order to describe den- 
sities, entropies and fluxes defined on the different levels 
of coarse graining. To this end we consider the level-2 
partitioning (cf. Fig. ||) of the cells two time steps after 
the initial coarse graining. Again primes denote quanti- 
ties evaluated after the first application of the mapping. 
Similarly, double primes are used for the values after 
two iterations. Note that also the transition probabilitcs 
Imi s m and r m and the widths and heights of the corre- 
sponding columns and strips receive primes now, due to 
their possible dependence on g m and w m (which evolve 
in time). 

In every cell m there are 3 2 strips labeled by a pair 
of symbols (p, q) with p,q £ {L, S, R}; p indicates the 
strip (L, S, R) of Fig. the point is located in after one 
time step (irrespective of the cell index), and q specifies 
its position in the strip (L, S, R) after the second time 
step. The strips (p, q) have height bs Piq 
and carry densities £>" n (p, q) and w'J n (p,q). Here, 

r f m , for q = R 

1 S m , for q = S (Bl) 
l m , for q = L 



is a shorthand notation for the height after one time step, 
and the subscript 




for q 
for q 
for q 



L 
S 
R 



(B2) 
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is used as a book-keeping device to indicate the cell where 
the points in a given strip come from. An analogous def- 
inition holds for the quantities s m . q without tildes. 

We call g m {p,q), g m (q), and g' m [g m (p,q)w m (p,q), 
g m (q)w m (q) , and g m w m ] the level-2, level- 1 and level- 
densities [kinetic-energy densities] after two time steps, 
respectively. By construction, the level-0 densities coin- 
cide with the cell densities considered in preceding sec- 
tions. The densities defined on different levels are related 
to each other since coarse-graining preserves their aver- 
age value in each cell: 



Here m p , q denotes the action of (B2) applied to m p : 
m Pt q — {m p ) q . In order to follow the time evolution 
of the entropies, we consider coarse-grained entropies de- 
fined with respect to densities of different resolution. The 

(i) 

level-i entropy S m is defined with respect to the level-i 
densities. Thus, Sm = S m — —abg rn \a.{g m w^ y ), while 
the level-one and level-two entropies at time r and 2t 
take the form 



Qm ^ ] s m Q g m (q) 



^2srn q , P g" n {p,q), (B3a) 



and 



i)' 



-ab 



v 



(p) ln \Q'm{p)w'm 7 (p) L (B5a) 



X] g m l9 em(9) tfl m(9) 
9 



<Sl 2) = -ab V sL J mq , P g" n {p, q) In g m (p, q)w" n 7 (p, q) 



/ j a m,q°rn q ,\ 



(B3b) 



(B5b) 



Based on the time evolution of the fields, relations be- 
tween densities at different instants of time can be calcu- 
lated. Taking into account the respective action (|2l]) and 
(E5j) of the mapping on the densities g and w, we obtain 



respectively. 



0m q) 



"m q ,q 
^rn^q 



S m q ,q Sm qtP ,p 



S m,q S m q -P 
W m (P> q) = w 'm q (P)(l + T Qrn') 

= w m (1 + Tq m )(l + rq m '). 



Qn 



(B4a) 



2. Entropy balance on level- 1 strips 

In order to obtain the entropy balance on level- 1 strips, 
we relate S$" to s£ ] ' ■ Making use of @, @, (§3) 
and ( |B5b[) , the level-2 entropy after two time steps S m 



(B4b) of cell to can be worked out as 
Sm' = -ab^s m ^ q s m ^ p g' m (p,q)ln[g' m (p,q)w' m ^(p,q)} 



- ab ^Z s ' mq , q Q mq ^-^r 



lb S 'm q ,q s m 9 , p ,pQm g , p H.Qm q (p) w 'm q (p) ~* i 1 + T <7»/) 



/\-7i 



(B6) 



where g' mq = J^ p s mq:P g mq (p)- 

■ j (2)" 

By using (B5a), Sm takes the form 



c(2)" 



-ab 



9 



s 'm q , q B mq l n + I] s 'm 9 , q S£ g ' + ah 7 



/ y "m q ,q^m q 
9 



e' m Jn(l + Tq m '), 



(B7) 



™,9 g 

which, after carrying out the summation over q, leads to the level-one generalization of the the entropy flux 



A o(l)' = o(2)" _ o(l)' _ _ h I V 5 



0) In Q m -\(p)w m -lljP)-= 



r m -i 



abl 



m+l ' 



l r 



a6 (r;„ + l' m ) s m , P Qm(P) ln \Q m ip) w 'm 1 {p)) + ab 



9 



'7 2^ S m„?em, ln (l + r <7m')- ( B 8) 
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The entropy balance for the level- 1 strips reads then 
AS£>' = S&>" - = XS^' + A e S^' (B9) 
with the level-one irreversible entropy production 

AiS&y = s$y' - s$" . (bio) 



Again it is related to the loss of information on the micro- 
scopic state of the system when applying a coarse-grained 
description. 

If the densities g m are uniform in the respective cells takes the form 
m q , i.e., they do not depend on the partitioning label 
p, the above relations coincide with Eqs. (|37]) and J45| ) 
[in which g m is eliminated via (p2|)j. Typically, however, 
there is a non- vanishing difference 



This flux is the two-time-step generalization of the en- 
tropy flux A e S m . 

In order to establish the entropy balance, the change of 
the coarse-grained entropy A^sffl is considered (with- 
out coarse graining after the first step) 

A(2)c(0) = o(0)" _ o(0) = 0(0)" _ c(2)" , c(2)" _ c(0) 
— ^ m *~^m ' ^ m 

= 55-^"+A< 2) ^». (B13) 



Thus, the two-step irreversible entropy change A^ Sm) 



A^fiW = A(2) ^ } - A< 2) S<? = Sir - S%". (B14) 



(Bll) 



Again it is identified as the loss of information caused by 
coarse graining. 

Note that this rate of entropy production can be ex- 
pressed as 



which only disappears in the macroscopic limit, as shown 
below after the discussion of the entropy balance for 
coarse-graining after every second time step. 



a(2)q(0) _ a.gW +A S {1} ' 



(B15) 



where AiS$ = S$' 



S { ri +1) ' (j = 0,1). When coarse 



graining is applied after each time step, the entropy pro- 
duction is AiSm +AiSin > . Thus, SAiSm) also amounts 
to the difference between the irreversible entropy produc- 
tion of the cases, where coarse graining is applied after 
In order to denote temporal changes taken with a time each and after ever y second time ste P> respectively. 



3. Entropy balance for coarse-graining after every 
second time step 



lag 2t, we assign the superscript ^ to A. A direct con- 
sequence of (B8) is that the entropy flux A^ Sm after 
two time steps is 

A(2)c(0) = _c(0) , c(2)" 

e ^ m — 1 

= -S (0) + ~S W +S {2) " 



4. Evaluation of SAiS., 



(i)' 



We recall that Q m {q) and w m (q) denote the fields after 
two time steps, when coarse graining is applied on the 
level-one partition. The level-one entropy after two time 
A e S$+A e Sgy. (B12) steps is 



(B16) 



S^" = ~ab ^2 Kn, q d m (q) In [Q m {q)w' m {q) 
i 

s' m r 

= -ab s ' mq , q Q mq ln -^T 1 + C, 9 ~ab Q mq In (g mq w m ^{l + rq m ) 

q ".8 q 

Using the form (|o|) of the coarse-grained entropy, we obtain 

s' 

S £)' = ~ab Y S m q , q Q mg ln -JT 1 + Y S 'm g ,q S m q + ab 7 S 'm t ,qQm, H 1 + T Qrn) 



(B17) 



where Sm\ is the coarse-grained entropy of cell m q evaluated after one time step. This form of level- 1 entropy can 
straightforwardly be compared with the level-2 entropy (B7), leading to 

A. <?«'- <?(!)" _ 9(2)" - a' rQ(°)' _ 9W1 = s A-9(°)+r ,A-9 (0) +7 l1 A-S' (0) 

L -^i°m ~ J m °m ~ / , °m q ,q\Pm q °m q J — °m<- i -l J m 1 ' m 1 2 m — 1 m+l 

= AiS^ — r m [AiS^ — AiS m _ 1 ] — l m [AiS$ — AiSU +1 ) + (r m _ 1 — r m )AiS < n J_ 1 + (Z m+1 — l m )AiS m+1 , (B18) 



1G 



where we have used that s m — 1 — r m — l m . Express- 
ing r m and l m in terms of driving forces and transport 
parameters (|3(]) and keeping only leading order terms in 
a, one obtains 



5A.S = rDdKAiS 1 ^) 
- r(d x v)A t S^ - 



0(r 2 ,a). 



(B19) 



After division by or, the right-hand side of (JBlSj) 
contains the rate of irreversible entropy production 
AiS^/ar and its spatial derivatives, because AiS^/ar 
has a finite macroscopic limit. Consequently, the differ- 
ence in the entropy production 8AiS/ (ar), vanishes in 
the macroscopic limit, when t — * 0. 

Since the temporal change of the level-0 coarse-grained 
entropy is (by definition) unaltered by changes of the pre- 
scription for coarse graining 



A(2)o(0) = o(0)" _ 

= <?(°)" - 

= as^' 



c(o) 

m 

C(0)' , c(0)' _ o(0) 

^AS«^ (B20) 



also the entropy flux, i.e., the difference of the change of 
the coarse-grained entropy and the irreversible entropy 
production, takes the same macroscopic limit, regardless 
of the procedure for coarse graining. 

Note that the calculations given in this Appendix can 
be generalized in a straightforward manner to account 
for averaging only after every n-th time steps on any fi- 
nite level k of the partitioning of the cells |39| . All these 
approaches only differ by terms that can be expressed as 
a product of r and spatial derivatives of the macroscopic 
rate of irreversible entropy production. For larger n and 
K, more and higher derivatives appear. 
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